Discovery of a GeV blazar shining through the Galactic plane 
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ABSTRACT 

The Fermi Large Area Telescope (LAT) discovered a new gamma-ray source near the Galactic 
plane, Fermi J0109-I-6134, when it flared brightly in 2010 February. The low Galactic latitude 
(6 = -1.2°) indicated that the source could be located within the Galaxy, which motivated rapid 
multi-wavelength follow-up including radio, optical, and X-ray observations. We report the results 
of analyzing all 19 months of LAT data for the source, and of X-ray observations with both Swift 
and the Chandra X-ray Observatory. We determined the source redshift, z = 0.783, using a 
Keck LRIS observation. Finally, we compiled a broadband spectral energy distribution (SED) 
from both historical and new observations contemporaneous with the 2010 February flare. The 
redshift, SED, optical line width. X-ray absorption, and multi-band variability indicate that this 
new GeV source is a blazar seen through the Galactic plane. Because several of the optical 
emission lines have equivalent width > 5 A, this blazar belongs in the flat-spectrum radio quasar 
category. 
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Introduction 



Variability in gamma-rays, as in other wave- 
bands, is an essential diagnostic for identifying 
sources and determining emission mechanisms. 
Gamma-ray variability has been established for 
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numerous extragalactic sources, all active galac- 
tic nuclei (AGN) and gamma-ray bursts. Within 
the Galaxy, pulsars exhibit periodic variability on 
short (milliseconds to seconds) time scales corre- 
sponding to their spin period (and are otherwise 
relatively steady), and binary systems exhibit pe- 
riodic variability on long (hours to days) time 
scales corresponding to their orbital period (Abdo 
et al. 2010a). 

Accreting binary systems have exhibited tran- 
sient outbursts during which their X-ray emission 
increases by a factor up to ~10^. One particular 
X-ray binary, Cygnus X-3 (a microquasar), has ex- 
hibited flaring activity in GeV gamma rays that 
is correlated with radio outbursts (Abdo ct al. 
2009b). Transient emission from Cygnus X-1 
has also been claimed by AGILE (Sabatini et al. 
2010). 

Flaring activity has not been established for 
Galactic GeV source classes other than binary 
systems. However, the Energetic Gamma Ray 
Experiment Telescope (EGRET) on board the 
Compton Gamma Ray Observatory did detect 
a bright gamma-ray transient near the Galactic 
plane, which had no radio-loud flat-spectrum ra- 
dio counterpart and had properties inconsistent 
with those of an AGN or isolated pulsar (Tavani 
et al. 1997). This source was never identified and 
has been interpreted as evidence for a new class 
of variable gamma-ray sources. 

The unexplained EGRET event has motivated 
an ongoing search for new gamma-ray transients 
near the Galactic plane, as well as for variability 
among established Galactic gamma-ray sources. 
Prior to the discovery reported here, three tran- 
sients near the Galactic plane (Hays & Cheung 
2009; Yasuda et al. 2009) have been detected with 
the Large Area Telescope (LAT) on board the 
Fermi Gamma-ray Space Telscope. None of them 
has yet been identified. 

On 2010 February l\ the LAT (Atwood ct al. 
2009) discovered a new, flaring source (Fermi 
J0109+6134) near the Galactic plane (Vanden- 
broucke & Hill 2010). The AGILE gamma- 
ray telescope confirmed the Fermi LAT detec- 
tion of the flare (Sabatini et al. 2010). We re- 
port an analysis of the full 19-month LAT data 
set for the source region. We also report new 
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Fig. 1. — Upper panel: Fermi LAT light curve, for 
photons with energy above 200 MeV. The source 
was not detected in a combined data set contain- 
ing all events between 2008 August 4 (MJD 54682) 
and 2009 May 30 (MJD 54971). This -10 month 
time interval is represented with a single upper 
limit. Lower panel: 2-10 keV X-ray light curve. 
There are seven Swift XRT observations (filled 
squares) and one Chandra observation (open cir- 
cle). X-ray fluxes are shown as measured, without 
accounting for absorption. Vertical error bars are 
±1 a uncertainty (statistical only). 

multi- wavelength observations, including an opti- 
cal spectrum from which we have determined the 
source redshift. Finally, we present a compila- 
tion of the broadband spectral energy distribution 
(SED) from both contemporaneous and historical 
data. The redshift, optical line width, SED shape, 
and multi-band variability indicate a blazar iden- 
tity for the source. 

2. Fermi LAT analysis 

We analyzed LAT data in a 15° radius around 
Fermi J0109-I-6134 for the full data set acquired to 
date (2008 August 4 through 2010 February 23). 
We performed an unbinncd likelihood flt to the 
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data using the gtlike tool, which is part of the 
Fermi Science Tools^. The statistical significance 
of the background-plus-source model relative to 
the background-only model is quantified with the 
"test statistic," TS = 2Alog(likelihood). The un- 
binned likelihood method is described in Mattox 
et al. (1996) and Abdo et al. (2010a). The back- 
ground model we used includes two diffuse com- 
ponents, one Galactic (model gll_iem_v02) and 
one isotropic (model isotropic_iem_v02), as well 
as nearby sources in the first Fermi LAT cata- 
log (IFGL; Abdo et al. 2010a). 

We fit Fermi J01G9+6134 with a power law en- 
ergy spectrum, $(£') = $o£^~'"^, with both the 
normalization $o and gamma-ray photon index 
free. The normalization of each of the following 
background components was also free in the fit: 
each point source in IFGL within a 6° radius of 
Fermi J0109-I-6134, the nearby bright and variable 
binary LS I -|-61°303, the Galactic diffuse model, 
and the isotropic diffuse model. Only those events 
belonging to the highest quality ("diffuse") class 
and with energies between 200 MeV and 300 GeV 
were included. Events with zenith angle greater 
than 105° were removed to avoid contamination 
by gamma rays produced from cosmic ray interac- 
tions in the atmosphere (Abdo et al. 2009a). We 
used the P6_V3_DIFFUSE instrument response func- 
tions. 

No statistically significant flux was detected 
above background at the location of Fermi J0109+6134 
in the 10-month LAT data set from 2008 August 
4 through 2009 May 30: The TS for the source 
hypothesis using all data in this time interval is 
3. The 95% confidence upper limit on the flux 
of photons above 200 MeV is 9.3 x 10^^ photons 

— 2 —1 

cm s . 

For the following 9-month interval, from 2009 
May 30 through 2010 February 23, Fermi J0109+6134 
was detected with a TS of 639, corresponding to 
a signiflcance of ^^25 a. The integral flux in this 
interval is (0.93 ± 0.05) x 10"'^ photons cm-^s'^ 
above 200 MeV, with = 2.59 ±0.06. Uncertain- 
ties are statistical only. According to the best-flt 
source model, 1577 photons with energy greater 
than 200 MeV are attributed to this source in 
the 9-month time interval. From this analysis. 
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the best-flt LAT position^ is {RA = 0l'^09'"58!4, 
Dec = -|-61°32'58"). The LAT 68% error radius is 
1.8', corresponding to a 95% error radius of 3.0'. 

No GeV source was previously known at this lo- 
cation. It was not detected in the EGR (Gasand- 
jian & Grcnier 2008), 3EG (Hartman ct al. 1999), 
or IFGL (Abdo et al. 2010a) catalog. The nearest 
source in IFGL is IFGL J0131. 2-1-6121, an unas- 
sociatcd source 2.55° away (Abdo et al. 2010a). 

The data from the 9-month interval were di- 
vided into monthly time bins, and a likelihood 
analysis was performed in each of them. The 
source was detected signiflcantly (TS>10) in each 
month. The gamma-ray light curve is shown 
in Figure 1. The source is clearly variable on a ~1 
month time scale, with two distinct flares peaking 
in 2009 August and 2010 February. The flux was 
at least 30 times greater at the peak of the flare 
than when it was not detected prior to the flare. 
The power law index was determined for each time 
bin and showed no signiflcant evidence of variation 
in the 9-month analysis. 

3. Counterparts 

We have identifled radio, optical, and X-ray 
counterparts of the new GeV gamma-ray source. 
Their positions are shown in Figure 2(a). 

3.1. Radio 

In archival radio data we identifled a likely 
flat-spectrum radio counterpart to the new GeV 
source: GT 0106-^613 (VCS2 J0109-h6133). This 
is the only radio source known within the 68% 
LAT error radius. This source was first dis- 
covered by Gregory and Taylor in their radio 
survey of the northern Galactic plane (Taylor 
& Gregory 1983; Gregory & Taylor 1986), per- 
formed with the NRAO 91-m telescope at 5 GHz. 
They determined that it had compact morphol- 
ogy, in comparison with other sources that fea- 
tured doublet or triplet morphology indicating 
one or two radio lobes in addition to a central 
object. A precise position for the radio source 
was later determined as part of the VCS2 very 
long baseline interferometry survey (Fomalont 
et al. 2003): {RA = 0l'^09™46!34439 ± 0^00013, 
Dec = +61°33'30'.'4557 ± 0'.'0003). 



^Wo use the J2000.0 epoch throughout. 
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Fig. 2. — (a) Gamma, X-ray, optical, and radio counterpart positions. The background image is a counts 
map from the Chandra observation. The Fermi LAT 68% error circle is indicated by the white arc. The 
Fermi LAT best fit position is outside the figure, at {I = 125.138°, b = —1.243°), in the direction of the white 
arrow. The two optical candidate counterparts are indicated by yellow (Paredes et al. counterpart) and red 
(IPHAS counterpart) circles. The two optical error radii are estimated at 0'.'25. The VCS2 radio position is 
indicated with a black '+,' and the best-fit Chandra position is indicated with a blue circle of radius O'.'S. 
Both the optical and Chandra position uncertainties are dominated by systematic frame offsets. The radio 
and X-ray positions confirm that the IPHAS source, rather than the Paredes et al. source, is the correct 
optical counterpart, (b) Reprocessed VLBA 8.6 GHz image from VCS2 (Fomalont et al. 2003) showing a jet 
extension of order 10 pc. The beam (2 mas x 0.7 mas at position angle = 82°) is indicated to the lower left 
and the contour levels begin at 1.5 mJy/bm and increase by factors of 2. 



The high-resolution VCS2 image resolved few- 
milli-arcsec (few-parsec, see redshift determina- 
tion below) jet extension (Fomalont et al. 2003). 
This image^ is shown in Figure 2(b). 

3.2. Optical 

Paredes et al. (1993) followed up the vari- 
able radio sources detected by Gregory and Tay- 
lor with optical observations. They detected a 
candidate optical counterpart of GT 0106+613 
with an / band magnitude of 19.4 and posi- 
tion {RA = 0l'^09™46!33, Dec = +61°33'29'.'l), 
1.4" offset from the VGS2 radio position. How- 
ever, without a spectrum they were unable to de- 
termine whether the source was Galactic or extra- 
galactic. 

In addition to the Paredes et al. observations 



*http: //astrogeo . org/vcs2/vcs2_cat .html 



performed in 1992, the source lies within both the 
IPHAS (Drew et al. 2005)^ and UVEX (Groot 
et al. 2009) optical surveys of the northern Galac- 
tic plane. The IPHAS images (obtained between 
2003 and 2006) show the Paredes et al. source, 
and they also show another optical source (IPHAS 
J010946.33-H613330.5) which matches the radio 
position better than the Paredes et al. candidate. 
This source has a position {RA = 0l''09'"46!33, 
Dec = -f61°33'30'.'5), O'.'l from the VCS2 radio 
position. The IPHAS and UVEX magnitudes are 
{i' = 19.7, r' - i' = 1.2) for the Paredes et al. 
counterpart, and (i' = 19.7, r' — i' = 1.5) for the 
IPHAS counterpart that is closer to the radio po- 
sition (Knigge et al. 2010). 



^http : //www . iphas . org 
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3.3. X-ray 



Following the gamma-ray flare, we obtained a 
series of seven observations of the source region 
with the Swift X-ray Telescope (XRT)^ between 
2010 February 3 and 2010 February 25, each with 
a duration between 5 and 10 ksec. We obtained a 
10 ksec Chandra ACIS-S3 (AXAF CCD Imaging 
Spectrometer) observation, performed'' on 2010 
February 21-22. 

In the initial Swift XRT observation (Cheung 
et al. 2010), two X-ray sources were detected in 
the fuU XRT field. One was outside the LAT 95% 
error circle. The other was a bright source posi- 
tionally coincident with the radio candidate, with 
(RA = 01^09™46?86, Dec = -H61°33'29'.'3) and a 
90% error radius of 4'.'4 according to the initial 
5 ksec observation. 

Figure 2(a) shows the X-ray counts map ob- 
tained with Chandra. In the Chandra field of 
view (larger than the zoom shown here), there are 
several sources detected. One matches the VCS2 
radio position well. The next closest is signifi- 
cantly dimmer and is 3.9' away from the radio po- 
sition, outside the LAT 95% error circle. The best- 
fit position of the Chandra source corresponding 
to the radio position is {RA = 0l''09'"46!338, 
Dec = -K61°33'30'.'42), with neghgible statistical 
uncertainty that is dominated by a frame offset 
uncertainty of ~0'.'5. There is no evidence for ex- 
tension or jet morphology in the Chandra image. 

The Galactic column density in this direction 
is estimated to be A'h = 5.5 x 10^^ cm~^ from 
to the Leiden/ Argentine/Bonn (LAB) 21 cm sur- 
vey (Kalberla et al. 2005). We fit the Chandra 
X-ray spectrum with a hard (Fx = l.Wt^is) ^^i- 
sorbed (A^h = 8.2^1:3 x 10^^ cm~2) power law, 
with = 40.3 for 61-3=58 degrees of freedom. 
It is not surprising that the column density esti- 
mated from the X-ray absorption is significantly 
larger than that estimated from the LAB radio 
survey data, particularly for this direction within 
10° of the Galactic plane. The X-ray estimate is 
likely more reliable than the LAB estimate due to 
uncertainty in the neutral hydrogen spin temper- 
ature, which is used to infer column density from 
the 21 cm data. 



^ Swift Target ID 31604. 

^Chandra Observation ID 11685, Sequence Number 11685. 
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Fig. 3. — Optical spectrum obtained with Keck. 
Note that extinction has not been corrected. 

The measured (absorbed) Chandra fiux is 
1.15 X 10~i2 erg cm-^s-i in the 2-10 keV 
band. The model (unabsorbed) fiux is 1.23 x 
10^^^ erg cm^^s^^ in the same band. The ab- 
solute flux uncertainty is ^ 25%. We also flt 
the spectrum including a thermal component (a 
kT=1.5 keV MEKAL spectrum) flxed at the newly 
determined (see below) redshift, z=0.783. Al- 
though some excess counts appear near the ex- 
pected Fe Ka line and near other lines corre- 
sponding to elements of intermediate mass, this 
thermal component was only detected at 1.4 a 
signiflcance. 

We analyzed archival Swift Burst Alert Tele- 
scope (BAT) data from 2004-2009 and found that 
the source was detected with a signiflcance of 
~4.5 a. We used this data to estimate the spec- 
trum in the 14-195 keV band, following the pro- 
cedure described in Ajello et al. (2008, 2009b). 
The resulting spectrum is accurate to the mCrab 
level (Ajello et al. 2009a). 

4. Optical spectrum and redshift 

We obtained an optical spectrum on 2010 
February 9 using the Low Resolution Imaging 
Spectrometer (LRIS) with the atmospheric dis- 
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persion corrector (Oke et al. 1995) at the Keck 
1 telescope. A l"-long slit was placed at a posi- 
tion angle of 12.5 degrees east of north in order to 
include both optical sources in a single pointing. 
The observations were performed at high airmass 
and as a consequence the image quality full- width 
half maximum (FWHM) was -1.2". The blue 
arm of the spectrograph used a 400 1 mm~^ grism 
blazed at 4000 A for a resolution of -7 A FWHM; 
the upgraded red side of the spectrograph used a 
600 1 mm~^ grating blazed at 7500 A, resulting in 
a spectral resolution of ~4.7 A FWHM. 

Data reduction was performed with the IRAF 
package using standard techniques: an optimal 
extraction algorithm to maximize the signal-to- 
noise ratio (SNR); wavelength calibration from sky 
lines; spcctrophotometric calibration and telluric 
subtraction from observations of G191-B2B; and 
visual inspection to clean cosmic rays. Multiple 
exposures were combined into a single spectrum, 
weighted by the SNR. 

The Paredes et al. candidate showed Balmer 
absorption at rest, indicating that it is a field star. 
The IPHAS counterpart showed a red continuum- 
dominated spectrum, with broad UV emission 
lines consistent with a redshifted blazar subject 
to Galactic extinction. 

The spectrum of the latter source is shown in 
Figure 3. Five emission lines were identified: C HI, 
Mg II, Ne V, H(5, and H7. From the spectrum we 
determine the redshift of this source to be 0.783. 
The equivalent width of each of the Mg II, Ne V, 
Hi5, and H7 lines exceeds 5 A. The signal-to-noise 
ratio of the C HI line is too low to estimate its 
width. Because the strongest emission line width 
was greater than 5 A, we classify this blazar as be- 
longing to the flat-spectrum radio quasar (FSRQ) 
category rather than the BL Lacertae (BL Lac) 
category, following the scheme of the Candidate 
Gamma-Ray Blazar Survey (Healey et al. 2008) 
and the first Fermi LAT AGN catalog (Abdo et al. 
2010b). 

A fit to the continuum in the Keck spectrum 
(between 380 nm and 800 nm) gives vF^ = 
8.3t^;^ X 10^13 gi-g at 5.45 x 10^'' Hz 

(550 nm) and a — 1 = 0.i1q4, where vFy cx 
j^-(q-i) g^j^^ ^Yy^ optical photon index is Topt = 
a + 1. Galactic extinction has been removed using 
the maps of Schlegel et al. (1998). The uncertainty 



on the normalization and index include 10% uncer- 
tainty in the extinction maps and 30% uncertainty 
in the absolute spectrophotometry due to slit loss 
and nearby source confusion. 

5. Radio, optical, and X-ray variability 

Taylor & Gregory (1983) determined that this 
radio source exhibited long-term (> 1 yr) and 
possibly short-term (< 1 yr) variability. Only 
a few percent of the sources in their northern 
Galactic plane survey exhibited variability, and 
they requested additional observations to deter- 
mine which if any were Galactic. Duric & Gre- 
gory (1988) established that the source exhibited 
short-term as well as long-term radio variability. 

In addition to historical (pre-2010) radio flux 
densities measured at several frequencies between 
0.365 and 22 GHz, a new radio spectrum, con- 
temporaneous with the 2010 February flare, was 
measured with the Effelsberg 100-m telescope on 
2010 February 4. The flux density was measured 
at several frequencies between 2.6 and 32 GHz. 
At 8 (22) GHz, the flux density was greater than 
archival values by '^20% (—50%) (Fuhrmann et al. 
2010), indicating radio flaring behavior contempo- 
raneous with the gamma-ray flare. 

There is marginal evidence for optical variabil- 
ity over the several years of IPHAS frames. The 
absence of the IPHAS source in the Paredes et al. 
observations (when the Paredes et al. counterpart 
was detected clearly) provides stronger evidence 
of optical variability. Both candidate counterparts 
were imaged using a 120 s exposure of the AC AM 
camera on the William Herschel Telescope (WHT) 
on 2010 February 5 to check the brightness at the 
time of the observed gamma-ray flare. The Pare- 
des et al. candidate had r' = 20.9 and the IPHAS 
candidate had r' = 21.2, each the same as deter- 
mined from the 2003-2006 IPHAS images (Knigge 
et al. 2010). 

Figure 1 shows the X-ray light curve compiled 
from the Swift and Chandra observations. In the 
2010 February 3-25 time interval, the X-ray flux 
shows no significant variability. 
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Fig. 4. — Spectral energy distribution. Historical data are shown in grey. Simultaneous data acquired during 
the February 2010 flare are shown in black. Extinction has been removed from both the optical and the 
X-ray spectra. The Fermi LAT spectrum was determined with a likelihood analysis analogous to the one 
described in the text, performed independently in each energy bin. For energy bins with no statistically 
significant flux detected {TS < 9), a 95% upper limit was calculated and is indicated with arrows. Historical 
radio data are from Douglas et al. (1996), White & Becker (1992), Condon et al. (1998), Becker et al. (1991), 
and Petrov et al. (2007). Dashed lines indicate a third-order polynomial fit to each bump. 



6. Spectral energy distribution and inter- 
pretation 

We combined both historical and contempora- 
neous radio, optical, X-ray, and gamma spectra in 
the SED shown in Figure 4. The two-bump shape 
is consistent with typical blazar SED's where the 
low-energy bump is interpreted as synchrotron 
emission and the high-energy bump is interpreted 
as inverse Compton emission. 



For reference, we fit a parameterization (third 
degree polynomial) to each bump, following the 
method of Abdo et al. (2009c). This fit gives a syn- 
chrotron peak at ~10^'^ Hz and a Compton peak at 
~10^^ Hz, making this a Low Synchrotron Peaked 
(LSP) blazar following the classification scheme of 
Abdo et al. (2009c). 

The redshift, SED shape, optical line width, 
and multi-band variability of this source indicate 
that it is a background blazar seen shining through 
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the Galactic plane. It provides an example of ex- 
tragalactic GcV source identification and charac- 
terization in the difficult region very close to the 
Galactic plane. Blazar catalogs generally avoid 
the region within ±10°of the Galactic plane due 
to the significant extinction and source confusion 
in this region. 

X-ray measurements of background blazars can 
improve our understanding of the Galactic inter- 
stellar medium. With improved statistics, the ab- 
sorption fit to the X-ray spectrum may constrain 
both the column density and composition in the 
direction of this source, independent of the LAB 
21 cm survey. Blazars are believed to have negli- 
gible intrinsic absorption (Perlman et al. 2005). 

This source provides the first example of suc- 
cessfully determining the identity of a transient 
GeV source discovered near the Galactic plane by 
following up with collaborative multi-observatory, 
multi-wavelength observations. The counterparts 
have been unambiguously identified thanks to the 
precise optical and X-ray positions. Monitoring 
of this blazar is ongoing with the Fermi LAT and 
other facilities. The three other Galactic plane 
transients detected with the Fermi LAT in its first 
19 months of operation remain unidentified. Ob- 
servation and analysis are ongoing to identify these 
transients and to detect new GeV transients near 
the Galactic plane. With its large field of view and 
frequent full-sky scans, Fermi is very well suited to 
discover more GeV transients across the sky and 
in particular in the Galactic plane. Fermi obser- 
vations of these transients will likely lead to new 
understanding of known GeV source classes as well 
as the discovery of new classes of GeV emitters. 
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